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Abstract 

We construct a special asymptotic solution for the forced Boussi- 
nesq equation. The perturbation is small and oscillates with a slowly 
varied frequency. The slow passage through the resonance generates 
waves with the finite amplitude. This phenomenon is described in 
details. 



Introduction 

The Boussinesq equation is one of fundamental equations of the mathematical 
physics. It naturally appears when one investigates the propagation of long 
waves into nonlinear media. The Boussinesq equation is the basic model for 
the description of surface waves PP,0,[H|; cross-axis oscillations of the balks 
and rods. It describes wave phenomenons in plasma and optical fibers 
We consider the forced Boussinesq equation 

U tt - U xx + a(U x )l + bll xxxx = e 2 T. 

The external perturbation has the different meaning for various applications: 
periodic influences on rods or an external pumping in optical media. The 
characteristic feature is the periodic nature of the external perturbation. 
Here we consider the such type of the perturbation 

U u - U xx + a{U x ) 2 x + e>yU xxxx = e 2 f(ex) exp{iS(e 2 x, e 2 t)/e 2 } + c.c, (1) 

where e is a small positive parameter; a, 7 are constants. The amplitude f(y) 
of the driving force is bounded and rapidly vanished function, when \y\ — > 00. 
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We investigate the simplest case S = S(e 2 t) = (e 2 t) 2 /2. It allows us to avoid 
the complicated calculations and saves the essence of the phenomenon. 

Before the resonance oscillations are realized under the perturbation. In 
a neighborhood of the curve t — the driving force becomes resonant and 
it leads to the change of the solution behaviour. The solution contains two 
waves of the finite amplitude after the slow passage through the resonance. 
In contrast to known papers related to the slow passage through resonance 
jEj, [Zj, [Ej in this situation we obtain waves with the amplitude of 0(1). This 
large increase of the order of the solution follows from the resonance on the 
zero harmonic. 

The paper has the following structure. In the first section we briefly de- 
scribe the result of the article and explain the observed phenomenon. The 
second section presents results of numerical simulations. Then we realize an- 
alytical asymptotic constructions. The forced oscillating solution of order e 2 
is constructed in the third section. The forth section contains the construc- 
tions in a neighborhood of the resonant curve. The behaviour of the solution 
after the passage through the resonance are described in the fifth section. 
All constructed asymptotic expansions are matched in the manner of |2|. 

1 Main result 

The formal asymptotic solution for (P) is constructed in the certain domain. 
This domain is divided on several subdomains and covers the resonant line 
t = 0. Our special solution has a different representation in each subdomain. 
The forced oscillations 

U (x, t, e) ~ e 2 u 2 (ex, e 2 t) exp{i(et) 2 / 2} 

describe the solution behaviour when t < 0. In a neighborhood of the reso- 
nant curve t = the asymptotic solution is represented by 

U(x,t,e) ~ w (ex,et) 

The dynamics of the leading-order term wq is described in terms of the Fresnel 
integral. After the passage through the resonance the solution becomes 0(1) 

U(x,t,e) ~ Vq (et + ex, e 2 t) + (et — ex, e 2 t) . 

This postresonant domain corresponds to t > 0. The leading-order terms 
Vq,Vq of the solution are determined from the pair of the Hopf equations 
(j2Hl)- All these representations are matched. 
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2 Numerical simulations 



In this section we present results of numerical simulations. They completely 
justify our asymptotic formulas and were done under the following conditions: 

1 f{eH)^ 



e = 0.05 and the perturbation looks like 



cos 



1 + (ex) 4 V 2e 2 
shows the generation of two waves of the finite amplitude of 0(1) 



. The figure 




3 The first external expansion 

In this section we construct the asymptotic solution of (JJJ that corresponds 
to forced oscillations. This solution has the order e 2 and oscillates under the 
perturbation. 



3.1 Equations for coefficients 

The WKB solution has the form 

U(x,t,e) = e 2 u 2 (xi, x 2) h) + e A u A (xi 1 x 2 ,t 2 ) 

+ e 6 u e (x 1 ,x 2 ,t 2 ) exp{iS (t 2 ) / e 2 } + c.c, 



(2) 



where x r 



e m t for m = 1,2. 
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Substitute this representation into (0) and gather the terms with the same 
order with respect to e. It yields the series of the algebraic equations 

(S') 2 u 2 = -/, (3) 

(S') 2 u 4 = 2tS'd t2 u 2 + iS"u 2 - d 2 Xi u 2 , (4) 

(S') 2 U(i = iS"u A + 2iS'd t2 u 4 + d 2 2 u 2 - d Xl u 4 (5) 

Note that 5" = t 2 in our model case. It allows one to obtain the exact 
representation for u n . 



U 2 = -if ( 6 ) 



n 



f 

4 

ihf-imij-dtj 



« 4 = J T 4 xlJ , (7) 



/ 6 

'2 



3.2 The domain of validity for external expansion (J^J) 

One can see that equations (jSJ) - © are not solvable, when 5" = 0. Moreover 
the representation (J2J) loses the asymptotic property in a neighborhood of 
the curve 

l[S] = S' = t 2 = 0. (9) 

Using (JHJ) - (JHJ) we obtain the representation of the solution in a neighborhood 
of the resonant curve 



U(x,t,e) = 

+ e 



{ 4 )\ exp V 2 ) + c ' c - (10) 

Representation (j2J) saves the asymptotic property when 



e max |{7 n +i| = o max \U n 

X\,X2M \X\,X2-t2 



It yields the domain of validity 

-t 2 > e. 
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4 The internal expansion 

In this section we construct the asymptotic expansion of the solution in a 
neighborhood of resonant curve (|9|). Here we use new scaled variables x\ and 
t\ and construct the solution of the form 

U(x,t,e) = w Q (x 1 ,t 1 ) + swi(xi,ti) + c.c. (11) 

4.1 Equations for coefficients 

Let us to substitute (fTT|) into (JTJ) and gather the terms of the same order 
with respect to e. 

The leading-order term satisfies to 

dlw - d 2 Xi w = f{ Xl ) exp{itl} + c.c. (12) 

The first-order correction term is determined from 

- d 2 xi w x = -2a [d Xl w d 2 Xi w Q + d xi w* d 2 xi w ] + c.c. (13) 

4.2 Solutions and asymptotics 

First of all we solve equation (fT2"j) . The equation for characteristics is 

fduV f du\ 2 . . 

UJ-fe) =0 ' (14) 

Introduce new characteristic variables 

y = t\ + x\, z = t\ — x\. 
Then we obtain the equation for v(y, z) = Wo(xi,ti) 

The general solution of (jT5j) is represented by 

Av(y,z) = V 1 (y)+V 2 (z)+j Z j" f {^^j {^(^2^)} dad ^ +c ^ 

(16) 

where Vi(?/) and V^z) are arbitrary continuously different iable functions. 
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To determine these functions V\(y) and V^z) we match this leading-order 
solution with expansion (j2J). Let us to return to original variables 

iw (ti + x 1 ,t 1 -x 1 ) = W+{t 1 + x 1 ) + Wv{t 1 -x 1 ) (17) 
1 f' 1 ^ 1 f tl+Xl , (a- p\ (ir^ + /3\ 2 \ J Ja 

and construct the asymptotics of (|T7|) . when t\ — > —00. The desired asymp- 
totics is obtained by integration by parts 



w 



ex P { y } J2 w ^ 2k + Wo + (*i + *i) + w-(T(*i - zi), (is) 
fc=i 



where w ,2 = -fM, w ,4 = 3z/ + d^/ 

Matching asymptotics (|18|) and ()10|) allows one to determine terms outside 
the integral in (fT7|). It yields W +(ti + ar x ) = W ~(*i - a?i) = 0. 

The solution of ()13|) has the form 

wi = j / ^(aj^dad/S + W^^i + arO + Wf^i-xi), (19) 

where 

tf(a, /3) = -2a (d^o^o + d^u^Wo) • (20) 
After substituting (JTHjl into (f2HJ) we obtain 



1 

Wl = - 



' "+/3 ^4 y 2 



dadf3 + + Wf + c.c. (21) 



Taking the integral by parts gives 



= + o(tr 8 )l ex P {^} + ( 4^ + o(t r 4 



+W 1 + (*i +W 1 -(ti -xi) +c.c. (22) 

In these formulas we keep the leading-order terms of expansions only. This 
rough approach is sufficient to obtain the domain of validity of (jllj) when 
ti — > —00 and match the solution with representation (fTUj) . 
Matching with JE} yields W?(ti +x t ) = Wf (*i - x x ) = 0. 
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Formulas (JTSJl and (|22|) allow us to determine the domain of validity for 
(fTTj) when t\ — > — oo. The condition 



emax |iu n+ i| = o I max |«j n | J , e — > 
is equivalent g<l. It means that expansion (jll)) is valid for t\ < 0. 
4.3 Asymptotics as £i — >• +oo 

In this subsection we construct asymptotics of w n as t 1 — ► +oo and determine 
the domain of validity for internal expansion (jlljl . 

Fist of all we investigate the leading-order term wq. As was shown above 
the coefficient Wq is represented by formula ()16|) . This formula can be written 
as follows 

w = - I I f ( - J exp <J i ( — ^— ] }> dad/3 



oo J — oo 



' ^ /( — i exp < i ( a + \ dadf3 



4 



oo ./ ii+xi 



/ / j ( ^ — — ) exp < % ( a + (3 ) \ dadf3 
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tj — aii ./ — oo 



J /-OO POO Q\ j / q, jj\ 2 

+ - / / ( J exp { i [ ) } dadp + c.c. 

After evident denotations we obtain 



w = C + Wq^x + xi) + w (t\ - Xi) H -j 1 exp \ f + c.c. (23) 



R{xiM) f.ti 

^ 6XP V2 

The last term of (j2SJ) decreases when \ti\ — >■ oo. This representation is 
obtained after single integration by parts. 

To obtain the asymptotics of W\ when t\ — > +oo let us to analyze right- 
hand side (|2()jl . It contains several character terms 



fi = g + (t 1 +x 1 )+g~(t 1 -x 1 )+g(t 1 + x 1 ,t 1 -x 1 )+ 1 J,' 1 exp{^/2} 
+ ^^-■Ijffl + ^.pi,,!)^, (24) 



t 2 
' 1 
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where 



g^h ± Xl ) = -2a(d xi w±d 2 xi w± + d xi w^ xi wfj, 

gfa + xuh-xi) = -2a(d Xl w^ d^WQ + 8 Xi Wq d 2 Xi w^ 
+ d xl w* + d 2 Xi WQ + d Xl w*~d 2 Xi w^ , 



Fi(x!,ti) = -2a{d Xl R d 2 Xl w^ + d 2 Xi w 

+ d 2 Xl R d Xl w£d Xl WQ + d Xl w* + d Xl WQ~ 

F 2 (x 1 ,t 1 ) = -2ad Xl R* [d 2 xi w+ + d 2 xl Wo 

F 3 {x 1 ,t 1 ) = -2ad Xl Rd 2 Xi R, 
F^xuh) = -2ad Xl R*d 2 xl R. 
After integration we obtain the following asymptotics 



Fi 

wi = (ti - xi)w+ + (ti + xi)w^ + wf + -j exp{itj/2} 
+5 exp{-it 2 /2} + ^ exp{i*i} + 5 + c.c. 



Here we denote 



w 



y rti+xi 

t = t / g + (y)dy, 

^ J -co 

ti—xi 



■y rti—xi 

w i = j / g~(y)dy, 

^ J —CO 
^ J —CO J —CO 



%^exp { ^/2} 



tl— xi /•ti+iri 



x exp 



i f y + z 



y+z y-z 
2 ' 2 

-2 



y + 2 



dydz 
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^^exp{- it ;/2} 



tl—xi pt±-\-xi 



oo J — oo 



y + g y - g 

2 ' 2 



X (>xp ~2~ (^2^ ' r ( ) ^' /: - 



(-1 



OO «/ — oo 



y + z y-z 



x exp < % 



y + z 



2 ' 2 

■i 



y + z 



dydz, 



t? "4y_ 00 M 2 ' 2 J I 2 



dydz. 



To determine the domain of validity for (jTJ) we use the inequality 



emax |w n +i| = o max |w n] i , 

X\,t\ \ Xi,ti 



6^0. 



It yields 



5 The second external expansion 

In this section we construct the second external expansion that describes the 
behaviour of the solution after the slow passage through the resonance. 

5.1 Equations deriving 

We construct the solution of the form 

U(x,t,e) = C + v£(tx + xi,t 2 ) +Uq(*i —Xi,t 2 ) 



+ e 



+ e< 



V^(tl + Xi,t 2 ) + V-l (ti - Xi,h) + Vi(ti + Xi,ti - X 1 ,t 2 ) 

v 2 ,s(x 1 ,t 2 )exp{iS/e 2 } +v£(ti + Xi,t 2 ) 



+ v 2 (tj - xi,t 2 ) + d 2 (*i + aci, *i - aci, ^) 



+ c.c. 



(26) 



The coefficients of (J2*H|) are determined recurrently. The function v 2t s is 
determined from the traditional algebraic equation 



-(S')Ws = f- 



(27) 
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On the next step we obtain the differential equation 

2d t2C v+ + 2ad c v+d 2 v+ + 2d t2V Vo + 2ad v v^d 2 ^ + Ad 2 ^ = 0, (28) 

where £ = t\ + x% and rj — t\ — x\. Consider this equation as equation for 
V\. It easy to see there is a special structure of the right-hand side of the 
equation 

$,«i = / + (C,fc) + / + fo,fc) 

To avoid the secular increase of V\ we obtain that vi = 0. It allows us to 
derive equations for the leading-order terms 



dl c v+ + ad c v+d 2 v£ = 0, 
d tr, v o + ad v v a 9 n v o = °- 

By substitution V + = d^v^ and V~ = 8 v Vq we get the pair of the Hopf 
equations 

d t2 V + + aV + d c V + = 0, 

d t2 V~ + aV~d v V- = 0. (29) 
Functions v f , ff, v%, v 2 are determined in much the same way 

5.2 The domain of validity for ( 1261) 

The usual condition for validity 

e max \v n+ i\ = o I max |i> n | J , e — > 

yields 

t 2 > e. 

Concluding remarks 

The realized rough analytical and numerical calculations allow one to obtain 
the finite amplitude waves due to the slow passage through the resonance 
of the small driving force. As it was mentioned above this large increase of 
the amplitude takes place due to the resonance on the zero harmonic. It 
becomes clear that this approach is evaluable for the phase S(x, t) of more 
general type. The construction and matching of infinite asymptotic series are 
possible also. The such type representation of the solution for the original 
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problem opens the way for a justification of the constructed asymptotics. 
But this is a theme for our future investigations. 
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